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QS^ The anharmonic contribution to phonon lifetimes and its temperature dependence 

^ is calculated from first principles in C, Si, and Ge using third-order density- 

|— { functional perturbation theory. Good agreement with available experimental data 

is obtained. Different competing two-phonon decay channels are compared and 

correlated with the density of final states. 
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The anharmonic decay of phonons into modes of lower energy determines the 
homogeneous width of first-order Raman lines and influences a number of other 
q important effects in semiconductors, including the decay rate of optical phonon 

populations as probed by time-resolved experiments. 1 Most of the theoretical 
studies performed so far in this field rely on simplifying assumptions on the 
0*\ possible decay mechanisms, 2 and/or on phenomenological models for describing 

harmonic and anharmonic interactions. 3 As discussed in Ref. 4, the ability of these 
schemes to account for experimental findings is however rather questionable. In 
i recent years, first-principles techniques based on density-functional theory have 

shown to be a very accurate and computationally viable tool for predicting various 
q vibrational properties of semiconductors. The use of these techniques to study the 

O anharmonic decay of phonons would therefore be of great interest. 5 Here we tackle 

this problem using third-order density-functional perturbation theory. 6-9 Our 
results are in very good agreement with some of the available experimental data, 
which however show a considerable spread. 

If only three-phonon processes are considered, energy and momentum 
conservation dictates that the LTO phonon decays into a pair of phonons with 
opposite momenta, whose frequencies sum up to the frequency of the decaying 
mode. The inverse lifetime of the LTO mode reads: 4 
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where n are the thermal occupation numbers, the j's indicate the phonon branches 
(j = 1 6), and V are anharmonic coupling coefficients, as defined e.g. in Ref. 4. 

Phonon inverse lifetimes, Eq. (1), have been obtained using phonon dispersions 
from Refs. 7 and 10, and anharmonic coupling coefficients calculated following Ref. 
9. The calculation of anharmonic couplings makes use of the so called In + 1 
theorem which states that the knowledge of the electronic charge-density response 
of a system up to order n in the strength of an external perturbation is sufficient to 
determine the energy derivative with respect to the strength up to order In + 1. 11 
For n = 1, this theorem ensures that third-order anharmonic couplings can be 
calculated from the linear response of the electron density distribution to lattice 
distortions, hence from the same ingredients which enter standard lattice-dynamical 
calculations in the harmonic approximation. 7 Our calculations have been performed 
within density-functional theory in the local-density approximation and using the 
plane-wave pseudopotential method. All the technical details are the same as in 
Refs. 7 and 10. 

Table 1. Calculated full widths 
at half maximum (2r) of zone- 
center optical phonons. The 
corresponding low-temperature 
experimental values are shown 
for comparison. The last col- 
umns indicate the relative con- 
tribution to r of the different 
decay channels. 

In Table 1 we report our calculated linewidths of the LTO phonons in C, Si, 
and Ge. The spread in the corresponding experimental data is rather large, ranging 
from 1.2 to 2.9 cm" 1 for C, 12 from 1.24 to 2.1 cm" 1 for Si, 13 and from 0.75 to 1.4 
cm -1 for Ge. 14 Among these, in Table 1 we display the values which are closer to 
our calculations, resulting in a very good agreement. 
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Figure 1. Temperature dependence of the full width at half maximum of the LTO 
phonon line in C, Si and Ge. Solid lines are the result of the present calculation; 
squares represent experimental data from Ref. 4. 



The temperature dependence of the phonon linewidth is given by Eq. (1) 
through the occupation numbers. The results are presented in Fig. 1 together with 
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the available experimental data. The agreement is very good up to temperatures 
T^>%ujlto I ks- Above these temperatures higher-order anharmonic terms are likely 
to account for the discrepancies. 

By restricting the sums over the j's to selected final states in Eq. (1), one can 
identify which are the relevant decay processes. The results of this analysis are also 
reported in Table 1, where the decay channels are indicated by 'TA' (j = 1,2) and 
'LA' (j = 3). The decay into one optical and one acoustic phonon is kinematically 
forbidden in all the three cases examined. For Si and Ge, the dominant channel 
involves one LA and and one TA mode as final states, consistently with the results 
of Ref. 5. The 'Klemens channel', i.e. the decay of the LTO mode into two acoustic 
phonons belonging to a same branch and with opposite momenta, turns out to 
give a very small contribution in these materials. The situation is quite different 
in diamond where the TA+TA channel becomes dominant at the expenses of the 
LA+TA channel. In order to obtain more detailed information on the relevant 
decay mechanisms, we decompose Eq. (1) into contributions, 7(0;), from different 
regions of the frequency spectrum. The function 7(0;) is defined by an equation 
similar to Eq. (1) where the sum over j\ and q is restricted to those values for 
which u^-^q) = u. This amounts to inserting 8 (a; — u^-^q)) under the sign of sum. 
Note that J ^(ui^dui = V. 
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Figure 2. Calculated phonon density of states, n(ui) (solid line), and spectral 
decomposition of the phonon linewidths, 7(w) (dashed line), for the three elemental 
semiconductors C, Si and Ge. 

In Fig. 2 we display 7(0;) and compare it with the one-phonon density of states 
(DOS), n(<jj). 15 By definition, 7(0;) is symmetric around wlto/2- At this frequency, 
7(0;) displays a peak corresponding to the Klemens decay mechanism. Lateral peaks 
may occur symmetrically with respect to ujlto /2 at frequencies uj\ and U2 such that 
n{ui) and n{uj2) are large and arise from a same region of the Brillouin zone. In 
diamond the Klemens peak is dominant because ujlto/2 falls in a region where the 
DOS is relatively large (between ^ta(L) an d ^ta(X))- The lateral peaks are less 
intense because the above conditions are not fulfilled in diamond. The opposite is 
true for Si and Ge, where n{ui,TO /2) is very small (and so is therefore the central 
peak). A careful analysis of the phonon dispersion relations reveals that the main 
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contribution to the lateral peaks comes from the decay into modes with j = 2 and 
j = 3 around the WL symmetry line, in agreement with the suggestions of Ref. 4. 

Our results show that third-order interactions provide an accurate description 
of the anharmonic decay of phonons in semiconductors up to far above room 
temperature. These interactions can be conveniently obtained from first principles 
by density functional perturbation theory with an effort comparable to that required 
for standard lattice-dynamical calculations in the harmonic approximations. This 
opens the way to the prediction of anharmonic lifetimes in systems — such as some 
bulk compound semiconductors (e.g. AlAs) or heterostructures — where they are 
not easily accessed by experiments. 
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